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Abstract —Line-of-sight (LoS) multiple-input multiple-output 
(MIMO) gives full spatial-multiplexing gain vrhen the antenna 
array geometry and orientation are designed based on the inter¬ 
terminal distance. These known design methodologies, that hold 
for antenna arrays with fixed orientation, do not provide full 
MIMO gains for arbitrary array orientations. In this paper, 
we study LoS MIMO channels with random array orientations 
when the number of transmit antennas used for signalling is 
2. We study the impact of common array geometries on error 
probability, and identify the code design parameter that describes 
the high signal-to-noise ratio (SNR) error performance of an 
arbitrary coding scheme. For planar receive arrays, the error 
rate is shown to decay only as fast as that of a rank 1 channel, 
and no better than SNR“® for a class of coding schemes 
that includes spatial multiplexing. We then show that for the 
tetrahedral receive array, which uses the smallest number of 
antennas among non-planar arrays, the error rate decays faster 
than that of rank 1 channels and is exponential in SNR for 
every coding scheme. Finally, we design a LoS MIMO system 
that guarantees a good error performance for all transmit/receive 
array orientations and over a range of inter-terminal distances. 

Index Terms —Antenna array, array geometry, coding scheme, 
llne-of-sight (LoS), multiple-input multiple-output (MIMO), 
probability of error. 


I. Introduction 

T he large swathes of raw spectrum available in the 
millimeter-wave frequency range are expected to provide 
an attractive solution to the high data-rate demands of the 
future 5G cellular networks HI. The small carrier wavelength 
of millimeter-wave frequencies allow for reduced spacing 
between the antenna elements when multiple antennas are 
used at the transmitter and receiver. This implies that multiple- 
input multiple-output (MIMO) spatial multiplexing gains can 
be obtained even in the presence of a strong line-of-sight (LoS) 
component when operating in such high frequencies IH. 

In LoS environments, the MIMO channel matrix H is a 
deterministic function of the positions of the transmitter and 
receiver and the geometry of the antenna arrays used at either 
terminals. If the positions of the communicating terminals are 
fixed and known apriori, the geometry of the antenna arrays 
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can be designed to optimize the performance of the commu¬ 
nication system. The LoS MIMO channel quality, in terms of 
capacity, multiplexing gain, coverage and channel eigenvalues, 
have been studied in |l2l-||9l as a function of the inter-terminal 
distance and the inter-antenna spacing of transmit and receive 
arrays, when the antennas are to be arranged in a rectangular, 
circular or a linear array. However, these design techniques 
assume that the position and the orientation of the antenna 
arrays are fixed, and the resulting criteria may be difficult to be 
satisfied if either of the communicating terminals is mobile or 
if the positions of the wireless terminals are not known a priori. 
Systems designed according to these known criteria degrade 
gracefully with variations in the geometric parameters, and 
may be adequate in certain scenarios where the changes in the 
orientation are limited, such as in a sectored communication 
cell where the variation of the base station orientation with 
respect to the direction of propagation is limited. However, 
these designs, which utilize two-dimensional antenna arrays, 
do not provide MIMO spatial multiplexing gains for arbitrary 
array orientations. 

In ITOl . the mutual information rates of a predominantly 
LoS channel with arbitrary antenna array orientations were 
studied using simulations and direct measurements in an in¬ 
door environment. The results show that the three-dimensional 
antenna arrays obtained by placing the antennas on the faces 
of a tetrahedron or a octahedron provide mutual information 
rates that are largely invariant to the rotation of antenna 
arrays in indoor LoS conditions. Previous studies of three- 
dimensional antenna arrays for wireless communications have 
mainly studied the capacity of the resulting MIMO system 
in a rich scattering environment. In HD a compact MIMO 
antenna was proposed which consists of 12 dipole antennas 
placed along the edges of a cube. A 24-port and a 36-port 
antenna were designed in na by placing antennas along 
the edges and faces of a cube. In ifTJj and iflTl . 6-port and 
16-port antennas were designed on a cube, respectively, and 
the performance of the MIMO system in terms of capacity 
and channel eigenvalues in a richly scattering environment 
were studied. The objective of HD-Hl has been to design a 
compact array by densely packing the antenna elements while 
exploiting the degrees of freedom available in an environment 
that provides abundant multipath components. 

To the best of our knowledge, there has been no prior 
theoretical study of LoS MIMO channels where the transmit 
or receive antenna array orientations are arbitrary, as may 
be experienced in wireless mobile communications. Lurther, 
all previous work have focussed on optimizing the mutual 
information rates of the MIMO channel. In order to achieve 
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the information theoretic limits, we need code design criteria 
based on an error performance analysis of the communication 
channel. In this paper, we consider LoS MIMO channels where 
the number of transmit antennas used for signalling is 2 and 
both the transmit and receive arrays have random orientations. 
We study the impact of the geometry of the antenna arrays on 
the system error performance and design a LoS MIMO system 
that guarantees a minimum channel quality and good eiTor 
performance for arbitrary transmit and receive orientations 
over a range of inter-terminal distances. 

We model the 2-transmit antenna rir-receive antenna LoS 
MIMO channel H using the upper triangular matrix R ob¬ 
tained from its QR-decomposition (Section |n|l. This allows us 
to derive bounds on pairwise eiTor probability and identify the 
code parameter that determines the high signal-to-noise ratio 
(SNR) error performance of arbitrary coding schemes in LoS 
MIMO channels. 

We show that for any planar, i.e., 2-dimensional, arrange¬ 
ment of receive antennas (such as linear, circular and rect¬ 
angular arrays), the rate of decay of error probability is 
similar to that of a rank 1 LoS MIMO channel whenever 
the receiver undergoes random rotations. Further, for some 
coding schemes, including spatial multiplexing IITsI - IITtI , the 
eiTor rate with any planar receive array decays no faster than 
SNR~^ even though the channel is purely LoS and experiences 
no fading (Section HIH i. 

We consider the smallest number of receive antennas nr = 4 
that can form a three-dimensional, i.e., non-planar, aiTange- 
ment, and derive bounds on error performance when they form 
a tetrahedral array. In this case, the eiTor probability decays 
faster than that of a rank 1 channel and is always exponential 
in SNR iiTespective of the coding scheme used (Section riV-Al l. 
We then design a LoS MIMO system with a good error per¬ 
formance for all transmit and receive array orientations over a 
range of inter-terminal distances by using a tetrahedral receive 
array and adaptively choosing two transmit antennas from a 
triangular/pentagonal array at the transmitter (Section lIV-Bl i. 
Finally, we present simulation results to support our theoretical 
claims (Section fVTl. 

Notation: Matrices and column vectors are denoted by bold 
upper-case and lower-case symbols respectively. The symbols 
AT, and HAUf denote the transpose, the conjugate- 
transpose and the Frobenius norm of a matrix A. The symbol 
II • II denotes the 2-norm of a vector. For a complex number 2 ;, 
arg( 2 :) and Re(z) denote its phase and real part, respectively. 
The expectation operator is denoted by lE(-). 


IT The 2 x LoS MIMO Channel 


We consider MIMO line-of-sight (LoS) transmission with 
rit = 2 antennas at the transmitter and Ur > 2 antennas at 
the receiver. Assuming that the large scale fading effects, 
such as path loss, are accounted for in the link budget, we 
take the magnitude of the complex channel gain between any 
transmit-receive antenna pair to be unity. If is the distance 
between the n}^ transmit and the m**' receive antennas, then 
the (to, component of channel matrix H G is ID 


= exp 


27rr„ 


( 1 ) 


where A is the carrier wavelength and i = s/—\. The re¬ 
sulting wireless channel is ypx = VSNRHx -f wrx, where 
yRx G C"’’ is the received vector, x G is the transmitted 
vector, wrx G C"'' is the circularly symmetric complex white 
Gaussian noise with unit variance per complex dimension, 
and SNR is the signal-to-noise ratio at each receive antenna. 
The power constraint at the transmitter is E (||x|p) < 1. We 
assume that the channel matrix H is known at the receiver 
but not at the transmitter. Let hi, h 2 G C"’’ denote the two 
columns of H, and H = QR be its QR decomposition where 
Q G has orthonormal columns, i.e., Q is a semi-unitary 

matrix, and 


llhill 

hlh2 

llhill 

0 



R = 


Let p denote the correlation between the two columns hi and 
h 2 of H, and 0^ be the phase of h|h 2 , i.e., 

|h|h2 


M = 


and 0fj, = arg (^h|h 2 ) . 


Jhill ||h2|, 

From ([T]), we have ||hi|| = ||h 2 || = yTv, and hence. 


R = \/tV 


1 e *®"p 


_o x/i - m". ■ 


( 2 ) 


Since Q is semi-unitary and wrx is a white Gaussian noise 
vector, y = Q Vrx is a sufficient statistic for x. Hence, in the 
rest of the paper we will consider the following equivalent 
channel 

y = V^Rx -f w, (3) 

where R is given in (|2]i, and w = Q^x is a two-dimensional 
circularly symmetric complex white Gaussian noise with zero 
mean and unit variance per complex dimension. 


A. Modelling the R matrix 

To analyze the error performance of arbitrary coding 
schemes in LoS MIMO channels, we model the phase as 
independent of p and uniformly distributed in [0,27r). Deriving 
the probability distribution of and p appears difficult, 
however, we provide an analytical motivation and numerical 
examples to support the validity of our model. 

We follow the notations from 0, 11 to describe the 
geometry of the transmit and receive antenna positions as 
illustrated in Fig. [T] We denote the inter-antenna distance at 
the transmitter by dt, and define the origin O of the three- 
dimensional reference coordinate system as the mid-point 
between the two transmit antennas. Define the z-axis of the 
coordinate system to be along the line connecting the two 
transmit antennas, i.e., the positions of the two transmit an¬ 
tennas are [O, 0, and [O, 0, —respectively. Choose 
the x-axis of the coordinate system such that the centroid O' 
of the receive antenna array lies on the x-z plane. Let O' be 
at a distance of R from O and at an angle [3 to the a:-axis 
i.e., at the point [i?cos/3, 0, i?sin/3]^. Consider an auxiliary 
coordinate system with O' as the origin and the three axes 
x',y',z' defined as follows: the x' axis is along the direction 


A 
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dt 



rxr, 


The distance 


between the n* transmit and receive 


th 


antennas satisfies Ed 


rm,n ~ i? + dm sin6»m coscjim + (-1)”— sin/3 + 

{dm sin 6>m sin0m)^ + (dm cosdm + (-l)”^COS/3)^ 

Therefore, the difference rm ,2 — fm.i is given by 

. ^ (dm cos dm + 4 cos/3)^ 

Cm,2 - Cm,l = d* Sin/ii H- — - 

ZH 

{dm COS dm - ^ COS ,5)^ 


= dt sin /3 + 


dtdm cos/3 cos d„ 


2i? 


F(/3)=hlh2= 

m=l 

'i2Trdt sin/3 


= exp 


A 


^exp 


m—1 


i2TTdtdm cos /3 cos d„ 


i?A 


Let /i(/3) = exp (i27rdt sin/3/A) and 

' i2TTdtdm cos (3 cos d„ 


ht) = XI 

m=l 


i?A 


We now upper bound the magnitude of the derivative of /i 
with respect to /3. The derivative of d/2/d/3 equals 

E "'' —i27rdtdm sin/3cosdm /i27rdtdm cos/3cosd„ 

-i?A-""P V -- 

m—1 ^ 


Note that |d/2/d/3| < 6, where b = 


YTm^i dr 


infinitesimal change A/3 in the value of , 


RX 


( 6 ) 

For an 


\f2{P + Ap)\-\h{P)\= /2{^) + ^A/3 


d/2 


Using the fact that 
we have 


kt + tu — kt 


|/2(/3 + A/3)|-|/2(/3)| 


< 


-|/2(/3)|. 

< Itul for any u,w £ C, 

d/2 


d/3 


|A/3| < 5|A/3|. 


Fig. 1. Illustration of the parameters used in the system model. 

OO', i.e., along the vector [cos /3, 0, sin /3] z' axis is on the 
x-z plane, and y' is parallel to y. Let {dm,dm,4>m) be the 
spherical coordinates of the wt receive antenna with respect 
to this auxiliary coordinate system, where dm is the radial 
distance, dm is the polar angle and (pm is the azimuthal angle. 


It follows immediately that |d|/2|/d,5| < b. Using the fact 
that p = \F{(i)\/nr = \f2{l3)\/nr, we have 


dp 

1 

d|/2| 

dp 

Tlr 

d/3 


< 


(7) 


Note that d^ = argU = arg/i + arg/ 2 , and hence, 
dd^/d/3 = d(arg/i)/d/3 + d(arg/2)/d/3. Now, arg/i = 
27rdt sin/3/A, and hence, d(arg/i)/d/3 = 27rdt cos/3/A. Us¬ 
ing O and the fact that the range of transmission R is much 
larger than dm, we have 


d(arg/i) 

d/3 


27rdtCOs/3 ^ 2Trdt I]m=i dr, 


A 


A 


Rrir 




dy 


d/3 


(4) 


Hence, we expect dd^/d,5 |d/i/d/3|, i.e., a small change in 

the value of /3, that causes a negligible change in y, changes 
the phase d^ by an entire cycle of 27r rad. This motivates the 
channel model where d^ is independent of y and uniformly 
distributed in the interval [0, 27r). 

Example 1. Consider a 2 x 2 LoS system operating in E-band 
at the frequency of 72 GHz over a distance i? = 10 m. Let the 
two receive antennas be positioned such that di = 0, d 2 = tt, 
pi = <p2 = d and di = d2 = dr/ 2 . Then, using (|5]l, we have 


Let U(/3) = h|h 2 denote the inner product between the two 
columns of H as a function of /3. Using ([TJ and (01, we obtain 


h{h 2 = 2 exp 
It follows that 
p. = cos 


f i2'Kdt sin /3 

V A 


cos 


Tvdtdr cos /3 

RX 


Trdtdr cos P 


and d^ = 


2'Kdt sin P 


( 8 ) 


(5) 


Then F(/3) = fi{P)f 2 {P), argU = arg/i -f arg/ 2 , and since 
I/ll = 1, we also have |F| = I/ 2 I. 

^The angle /3 is equal to the parameter 6t used in 


RX J X 

Suppose the antenna geometry is to be configured so that H 
is unitary, i.e., p = 0, under the assumption that /3 = 0. This 
can be achieved by choosing dt and dr so that 

dt dr cos P dtdr 1 
RX ^ 'RX ^ 2' 

This is the criterion for uniform linear arrays 
given in B-El. With A = 4.2 mm, the choice of 
dt = dr = ^/RX/2 = 0.145 m yields p = 0. With this 
choice of dt and dr, through direct computation using d^, we 
observe that as /3 undergoes a small variation in value from 
0 rad through 0.029 rad (1.66°), the corresponding value of 
p changes from 0 to 6.6 x 10“"^, while d^ ranges over the 
entire interval from 0 to 27r rad. ■ 
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Fig. 2. The joint probability density function of Example 


Fig. 3. The joint probability density function of Example 


Example 2. Continuing with the 2x2 system of Example [T] 
now assume that the transmit and receive arrays are affected by 
independent random rotations about their respective centroids. 
The random rotations are uniformly distributed over the space 
of all 3-dimensional rotations. The channel matrix H, and the 
parameters 0^ and p are now random variables. The joint 
probability density function obtained using Monte- 

Carlo methods is shown in Fig.|2] We computed /(6*^, p) over 
a rectangular grid of 625 points using 10^ randomly generated 
instances of H. For any fixed p, we observe that /(6>^,/i) is 
essentially constant across all values of 0^, implying that 0^ 
is uniformly distributed in [0, 27r) and is independent of p. ■ 

Example 3. Consider a 2 x 4 FoS MIMO system, with a 
rectangular array at the receiver, carrier frequency of 72 GHz, 
and inter-terminal distance of i? = 10 m. The receive antennas 
are placed at the vertices of a square whose edges are of 
length dr- We choose dt = dr = \J which yields 
the ideal channel (i.e., /i = 0) if the transmit and receive 
arrays are placed broadside to each other 13 . The joint 
probability density function /(0^, /i), obtained using Monte- 
Carlo methods, when the transmit and receive arrays undergo 
uniformly random rotations about their centroids is shown in 
Fig. [3] As in Example |2] the numerical result supports the 
validity of our channel model. ■ 


In the rest of the paper we model the 2 x Ur FoS channel 
using the 2x2 matrix (cf. ©) 


R = vTv 




0 y/1 - 


(9) 


where 0 is uniformly distributed in [0,27r) and 



771=1 


exp 


f i2T:dtdm cos j5 cos dr 




RX 


( 10 ) 


B. Coding schemes 

We analyse the error performance of any arbitrary coding 
scheme for two transmit antennas with a finite transmission 
duration. Fet T > 1 denote the transmission duration of a 
given communication scheme and C the finite set 


of all possible transmit codewords. The rows of the code¬ 
words X G correspond to the two transmit antennas and 
the columns to the T time slots. All codewords are equally 
likely to be transmitted and the optimal decoder, i.e., the 
maximum-likelihood (ME) decoder, is used at the receiver. 
We further assume that the communication scheme satisfies the 
average power constraint II^IIf — \‘^\ Our analysis 

holds for arbitrary codes including space-time block codes 
(STBCs) ifT^ . 

We now briefly recall two specific coding schemes which 
will be used in our simulations (in Section Iyj to illustrate 
our analytical results. Spatial multiplexing (SM) ca-Ei, 
which is also known as VBIAST in the literature, is a simple 
yet powerful scheme where independent information symbols 
are transmitted across different antennas and time slots. The 
codebook ^ C corresponding to SM occupies T = 1 

time slot, and is given by 





Si, S2 



where ^ is a complex constellation, such as QAM or PSK. 

The Golden code ED is an STBC for two transmit antennas 
occupying T = 2 time slots, and is given by 


J a(si-|-TS 3 ) a{s2 + TS4) 
\ ia{s2 + PS4) a{si + P.S3) 


Sl,.. 


, S 4 


G A 


where ^ is a QAM constellation, r = (1 -f •\/5)/2, p = l/r, 
a = 1 + ifj, and a = 1 + ir. Unlike SM, the Golden code 
spreads the information symbols across time and antennas. 

Both SM and Golden code have been well studied in the 
case of non line-of-sight MIMO fading channels. The SM 
scheme provides high data rate with low complexity encoding 
and decoding, while the Golden code provides high data rate, 
full-diversity as well as a large coding gain at the cost of 
higher decoding complexity in fading channels. 


C. Error probability analysis for a fixed p. 

We now analyse the error performance of a given arbitrary 
coding scheme for a fixed value of p. Fet ^ C be 
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any code and Xa,Xb G be two distinct codewords. Let 
AX = Xa — Xft be the pairwise codeword difference matrix. 
The pairwise error probability between X^ and X^ for a fixed 
/i and a given realization 0 = 6* is mi 


is the modified Bessel function of the first kind and zeroth 
order. For large x we have ll20l 

loix) = + O {x~^)) . (15) 

V^TTX 

Using (foi l and the first order approximation (fTSl l. we get the 
following approximate upper bound when /r > 0, 


where Q is the Gaussian tail function. Using the Chernoff 
exp (—x'^/2) 

bound Q{x) < -we have the upper bound 

PEP<iexp(^-^||RAX|||.). (11) 

Denoting the two rows of the matrix AX as Ax[ and Ax^, 
we obtain the following expression for the squared Euclidean 
distance between the codewords at the receiver, 

||RAX||^ = n,(||Axif + IIAxaf + 2 pRe(e*®Axl Axs)) 

= nr.(^||Axi|p + ||Ax 2|P + 2pcos0'|Ax|Ax2|^ 

( 12 ) 


where 0' = 9 + arg(Ax|Ax 2 ) mod 27r. 

1) Worst-case Error Probability over 9: For a given p, 
the value of 9 that minimizes the squared Euclidean distance 
||RAX|p at the receiver is 0* = tt + arg(Ax|Ax 2 ) since it 
leads to cos9' = —1 in (fTSll . Using the notation 

d(p, AX) = ||Axi|p + ||Ax 2|P - 2/r|Ax|Ax2|, (13) 


the worst-case squared Euclidean distance is 

min ||RAX|||, = Urdu, AX). 

SG [0,271-) 

Thus the worst-case PEP for a fixed p satisfies 

^ 1 f-nrSmd{n,AX)\ 

PEP*(m) < - exp f ^. ( 14 ) 

2) Average Error Probability over Q: Since 0 is uniformly 
distributed in [0,27r), so is 0' = 0 + arg(Ax| AX 2 ) mod 27r. 
Using (fTTT i and (fTSl l. the error probability averaged over 0, 
for a fixed fi, can be upper bounded as follows 


Ee (PEP) < Ee ( ^exp ( -r:A2||RAX||| 


SNR, 

T 


1 /-SNRn,(||Axif+ ||Ax2f) 


= -exp 


exp - 


SNRn^ 


2)rcosd'|Ax|Ax2| ) d9' 


= iexp 


-SNRn,(||Axi||2 + ||Ax2|p) 


/SNRrir I . t A A 
— 2 —At|AxjAx2|j 


where 


1 ^TT nZ-K 

Io{x) = — / exp (xcosd') d0'= ■— / exp(xcos0') d0' 
tt Jq 2 tt Jq 

p 2 -K 

/ exp (—xcosd') d0' 

Jo 


1 

2?! 


Ee(PEP)< - 

y 47mrSNR/r|Ax|Ax2| 

X exp AX)^ . (16) 

Since the exponential function falls more rapidly than 
SNR^^^^, the high SNR behaviour is dictated by d(/i, AX). 

In this section, we derived bounds on PEP for a fixed /i. In 
Sections flUl and flVl we analyze the effects of random rotations 
of the terminals on /i and error performance. 

III. Error performance oe planar receive arrays 

Assume that the receive antenna system is affected by a ran¬ 
dom three-dimensional rotation U G about its centroid 

O'. Let the rotation U be uniformly distributed on the set of 
all 3-dimensional rotations, i.e., the special orthogonal group 

SO 3 = {U G I UUT = I, det(U) = 1} . 

In Theorem [T1 we provide a lower bound on the average 
pairwise error probability over a LoS MIMO channel with 
planar receive array. To do so, we derive a lower bound on 
the probability that a random rotation U would lead to a ‘bad’ 
channel matrix with fi close to 1, i.e. p, > 1 —e for some small 
positive e. By analyzing the PEP for this class of bad channels, 
and letting e decay suitably with SNR, we arrive at a lower 
bound for the average PEP at high SNR. 

Theorem 1. Let the receive antenna array be any planar 
arrangement of Ur antennas, > 2, undergoing a uniformly 
distributed random rotation U about its centroid. At high SNR, 
for any transmit orientation /3, we have 

exp (_ "7 .c|AxIax.| 

E(PEP) > - ^ - 

2n,SNRy2^2|AxtAx2| (||AX||;^ + 

X exp d(l, AX)^ , (17) 

where c = 27rdtdm/i?A. 

Proof: Let {ea,, ej,, e^} be the standard basis in R.^. When 
the receive system undergoes no rotation, i.e., when U = I, let 
the position of the receive antenna relative to the centroid 
O' of the receive antenna system be drrJtm, where G R^ is 
a unit vector. Since the receive array is planar and the random 
rotation U is uniformly distributed, without loss of generality, 
we assume that the vectors ri,..., r„^ are in the linear span 
of Bx and e^. From Fig. [T] we see that 9m in (O is the angle 
between the orientation Ur^ of the mJ^ receiver and the unit 
vector V = [—sin/3, 0, cos/3]^ along z'-axis, i.e., cos 9m = 
rJ^'LP'v. Note that Ut has the same distribution as U, and 


) 
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V = U‘'’v is uniformly distributed on the unit sphere in 
The resulting random variable |ejv| is known to be uniformly 
distributed in the interval [0,1]. 

For a small positive number 5 > 0, consider the event 
8 : |ejvp > 1 — (5^. The probability of £ is 

P(f) = P (lejvl > x/1 - 52 ) 

for small values of <5. We will now derive an upper bound 
for the PEP for the case when £ is true. Using the fol¬ 
lowing inequalities, we first show that | cos 6 m\ ^5, for all 
TO = 1, , Ur, 

|cos6»„|^ = IrT^vp 

< |ej;vp -f lejvp (since S span(ea;, e^)) 

= l|vf-lejvp 

< 1 - (1 -( 5 ^) = 

Let Cm — COS/3/-RA 3.n(l Cmax — II13,x{ci, • ■ • , 

From (fTOb , we have 



Ur 

exp(zCmCOS0m) 

m—1 


We will now show that the value of p is close to 1 when £ is 
true. If Cm = 1 — exp(zcmcos^), then 

kmp = (1 - COs{CraCOs9m)f + sin^ [CmCOsOm) 

= 2 — 2 COs(Cm COS Om) 

cl^COS^{ 6 m) \ 

2 ) 


2-2 1 - 


where the approximation follows from the Taylor’s se¬ 
ries expansion of the cos{-) function and the fact that 
IcmCOS^ml < Cm <5 is Small. Now, 


M = — 

Tlr 


(1 Cm) 


1 

Ur 


Ur 


-E' 


^i--Ei 

n.. A—^ 


> 1 — 6 c„ 


Thus /r > 1 — (5cmax whenever £ is true. 

The pairwise eiTor probability for fixed p and 0 = 0 is 
Q(VSNR||RAX||y2y Since we need a lower bound on 
the probability of error, we use the following lower bound for 
the Gaussian tail function ll2Tll 

2 f x^\ 

Q(x) > , _, - , , exp (-I, for X > 0. 

Using x^ -f 4 < (x -f 2)2 for x > 0, we obtain a more relaxed 
bound 

Q(x) > exp . 

^ ^ - v^(x+i) ^ V 2; 


In our case x = a/SNR||RAX|||./ 2, and we use the exact 
value of X from (fTSl i for the exponent, and the following upper 
bound for the denominator 

X = yM||RAX||^ < ||R||^ IIAX||^ 

= \AySNR||AX||F. 


Thus, we have the following lower bound for a fixed p, and 

0 = 0 , 


PEP > 


exp(-^||RAX||2,) 

(VnrSNR||AX||F -f 1) ■ 


(18) 


Since the denominator is independent of the phase 0, we can 
use the same method as in Section ITI-C2I to obtain the average 
of the above lower bound over the uniformly distributed 
random variable 0 . Averaging (fTSl) over 0 and using the 
approximation to the Bessel function (fTSl) . we obtain 


Ee(PEP)> 


exp(-^d(p,AX)) _ 

n^SNRy27r2p|Ax|Ax2| (||AX||f + 7 ^^) 


Using the trivial upper bound p < 1 in the denominator, 
exp(-^d(p,AX)) 


Ee(PEP) > 


n^SNR^27r2|AxlAx2| (||AX||f + 


(19) 


Since d(p, AX) is a decreasing function of p, if £ is true, 
the numerator in the RHS of (fT9l l can be lower bounded 
by exp (—— ( 5 ci„ax 7 AX)). The expression ( fT9] l is 
a lower bound on the average PEP for a given p. We now 
derive a lower bound for the PEP when averaged over both p 
and 0 as follows 


E(PEP) = P(f)P (Xa ^ Xb|f) + P(£:")P (Xa ^ Xb\£^) 

> P(f)P(Xa^Xb|f) 

^ ^2 exp (-^d(l-fcmax,AX)) 

■ 2rr.SNR^27r2|AxlAx2| (||AX||^^ -f 

( 20 ) 

From the definition ( fTSl l of d(p, AX), we have 

d(l - JCmax, AX) = d(l, AX) + 2JCmax|AxlAx2|, 


Using the above relation and choosing <5 = SNR which is 
small for high SNR, we obtain 


E(PEP) > 


exp exp (-2^d(l, AX)) 

2n,SNRy2^2|AxlAx2| (||AX||f -f 


Using cos /3 < 1 in Cm = 2Trdtdm cos /3 / RX we obtain Cmax > 
maxm 2TTdtdm/RminX. This completes the proof. ■ 

We compare the lower bound (fTTI i on PEP for planar 
receive arrays undergoing random rotations, with the upper 
bound ( fTbl l for a channel with fixed p = 1 . The dominant 
term dictating the rate of decay of eiTor probability for both 
these channels is exp minAx d(l) AX)), where the 
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minimization is over all non-zero codewords difference ma¬ 
trices AX = Xq — Xf, of the code Note that /r = 1 
minimizes the performance metric d(/i, AX), and corresponds 
to the worst-case scenario where both H and R have rank 1. 
While planar receive arrays, such as the well-studied linear, 
rectangular and circular arrays, provide an array gain (an n^- 
fold increase in received SNR), their asymptotic coding gain 
minAxd(l, AX) provides no improvement over that of any 
rank 1 channel. 

Theorem[T]further implies that when miuAx d(l) AX) = 0, 
the error probability is no more exponential in SNR, but 
decays at the most as fast as SNR“^. Hence, although the 
channel is purely LoS and experiences no fading, the error 
performance with a planar arrangement of antennas can decay 
slowly, similar to a fading channel. 

The parameter d(l,AX) satisfies the following tight in¬ 
equality 

d(l, AX) = ||Axi|p -f ||Ax 2 |P - 2|Ax|Ax2| 

>||Axif+ ||Ax 2 f- 2 ||Axi|| ||Ax 2|1 

= (||Axl||-||Ax2||)^ (21) 

The second line follows from the Cauchy-Schwarz inequality 
which is tight if and only if Axi and Ax 2 are linearly 
dependent. Thus, d(l, AX) = 0 if and only if Axi and Ax 2 
are linearly dependent and ||Axi|| = ||Ax 2 ||, i.e., if and only if 
Axi = aAx 2 for some complex number a of unit magnitude. 
We use this observation in Example |4] below to show that the 
widely used spatial multiplexing coding scheme suffers from 
such a slowly decaying error probability with planar receive 
arrays. 

Example 4. Performance of Spatial Multiplexing with Planar 
Receive Array. The codeword difference matrices of the SM 
scheme are of the form 



where Asi, As 2 S A^ and AA = {x — y\x,y G A} is the 
set of pairwise differences of the complex constellation A. 
When Asi = As 2 the two rows of the codeword difference 
matrix AX are equal resulting in d(l, AX) = 0. Hence, for 
the SM scheme, miuAx d(l, AX) = 0, and from Theorem [T] 
the rate of decay of the average error probability will be no 
faster than SNR~^. Note that this result is valid for any number 
of antennas Ur used in any planar arrangement of the receive 
array. This theoretical result is validated by our simulations 
(see Fig. [TO] and Fig. fTsT i in Section |V] ■ 

IV. Error Performance oe Tetrahedral Receive 

ARRAY 

The smallest number of antennas that can form a non-planar 
arrangement is 4. In this section we consider the case where 
rir = 4: receive antennas are placed at the vertices of a regular 
tetrahedron, see Fig. @1 The inter-antenna distance dr is the 
same for any pair of receive antennas, and this is related to 
the distance dm of each antenna from the centroid O' of the 


1 



Fig. 4. The receive antennas are placed at the vertices of the 

tetrahedron. Also shown in the figure are the centroid O', the distances 
and ^4 of the antennas 3 and 4 from O', and the inter-antenna distance dp. 


receive array as dm = s/Wdr , m = 1,..., 4. Fet us define 
the deviation factor rj as in El, El as follows 


V = 


RX 

2dtdr cos/3 


( 22 ) 


In the case of a tetrahedral receiver, using (fTOl l and (l22l l. 



4 

H exp 

m—1 



The parameter 77 captures both the distance R and the transmit 
orientation /3, while the variables 61 ,. ■ ■ ,(^4 jointly determine 
the receive orientation U. In order to upper bound the error 
probability using (fT4l i. we need the maximum value of y, over 
all possible rj and U. Fet 


P*{ri) 


1 

max — 
UG503 4 


4 

H exp 

m—1 



(23) 


be the maximum channel correlation over all receive orien¬ 
tations as a function of rj. If one is aware of the range of 
values that R and (3 may assume, then one can upper bound 
the worst-case PEP using (fT4l) as 


PEP* < i exp 
= iexp 


Tlr SNR ^ ^ 

SNR d(maxp*(? 7 ), AX)^ . 


(24) 


A. An upper bound on y*(jf) 

In this sub-section we derive an upper bound on p* ( 77 ) for 
all 77 > 1. This result will allow us to show that the high SNR 
error performance of the tetrahedral array is better than any 
planar receive array when 77 > 1 and the receiver undergoes 
a uniformly random rotation. To derive this upper bound, we 
first show that when 77 > 1 , irrespective of the receive array 
orientation, the 4x2 channel matrix H contains at least one 
2 x 2 submatrix Hsub such that the correlation psuh between 
the two columns of Hgub is at the most cos {yx This 
latter problem is equivalent to finding the maximum distortion 
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Fig. 5. The tetrahedron arrangement illustrating the vertices the 

reference O' at the centroid of the tetrahedron, and the directions of a few 
of the unit vectors Vm and gm,£- 


when a unit vector in is quantized using a codebook Q 
consisting of 12 unit vectors that correspond to the 6 edges of 
the tetrahedron along with the polarities ±1. The computation 
of this maximum distortion is then simplified by showing that 
G is a group code ll22]| . 

We first introduce some notation to capture the geometrical 
properties of the tetrahedral array. Consider the tetrahedron 
shown in Fig. |3 with the centroid O'. Let G be 
the unit vector in the direction of the receive antenna 
with respect to the reference O'. Hence, the position vector 
of the receive antenna is dmYm- If one applies a 3- 
dimensional rotation U G R^^^ on the receive system about 
O', the position of the receive antenna is dmUr^,. It 
is straightforward to show that the polar angle 9m of the 
rotated receive antenna (cf. Fig. [TJ satisfies cos 9m = 
where the unit vector v = [— sin/3, 0, cos/3] Since U is an 
arbitrary rotation matrix, the set of all possible values assumed 
by the vector v = U‘'’v is the sphere consisting of all unit 
vectors in R^. From (fTOl i. the correlation p for a tetrahedral 
receiver is 



^exp 


f i2TTdtdm cos ,8 cos 


where cos 9m = ^nd v G captures the 

effect of the rotation undergone by the receive array. For any 
m ^ the unit vectors Vm and rg satisfy ||rm — r^jj = y^8/3. 
Let 




I’m - re 

\\rm - ri\\ 



-re) 


be the unit vector along — re, i.e., along the edge of the 
tetrahedron between the vertices m and £ (see Fig. |5]l. 

Let Hsub be the 2x2 submatrix of H formed using the 
and £*^ rows. Note that Hgub is the channel response 
seen through the receive antennas m and £. Using the fact that 


dm = de = 3/8 dr, the correlation between the columns of 

Hsub can be written as 


_ 1 
Msub — 2 


exp 


i2Trdtdm cos 13 


RX 


exp 


i2TTdtde cos (3 rjv 



^ i2TTdtdm cos j3{Ym - re)'' 

RX 

i2'Kdtdm\/W^ cos/3g^ 
EA 


cos 


(25) 


where the fourth equality follows from (l 22 T i and the last 
equality uses straightforward algebraic manipulations. Given 
an ‘orientation’ v, we intend to find the submatrix Hsub with 
the least correlation psuh- If ty > 1, we have 


— K 

2 r] 


TT 

< -. 

- 2 


l|v-g» 


Since cos is decreasing function in the interval [0,7r/2], 
from (l25l l. the problem of finding p^uh translates to finding 
the edge gm,i of the tetrahedron that has the largest inner 
product with v. 

We will now show that for any v G there exists a gm,e 
such that y^l/2 < < 1. Since 

= l|vf + ligm.rf - 2g;^_^v = 2(1 - 

this is equivalent to finding the maximum squared Euclidean 
error when the set of vectors G = {gm,i \ m y £} is used as 
a codebook for quantizing an arbitrary unit vector v in R^. 
The set G contains 12 vectors, corresponding to the 6 edges 
of the tetrahedron together with the polarity ±1. 

Proposition 1. For any v G there exist m,£ G {1, 2, 3,4}, 
m y £, such that g^ > \A72- 

Proof: With some abuse of notation we will denote the 
elements of G as gi,..., gi 2 . For each i = 1,..., 12, let 


A = {v G I gjv > gjv, for all j y i} 


(26) 


be the set of unit vectors that are closer to gi than any other 
gj G G- Since UiVi = S^, it is enough to show that 


min min gJv = 

I vGCi 

As we now show, the regions Vi,..., /D 12 are congruent to 
each other. Let FL be the symmetry group of the tetrahedron, 
i.e., the set of all orthogonal transformations on R^ that map 
the tetrahedron onto itself. It is known that the group FL 
is isomorphic to the symmetric group 184 of degree 4, and 
every element of FL is uniquely identified by its action on 
the set of vertices, which is isomorphic to the action of the 
corresponding element in ^4 on the set {1,2, 3,4}; see ||23]| . 
Since for any two given pairs (mi,£i) and {m 2 ,£ 2 ), with 
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Fig. 6. An illustration of the cones S and 7^i used in the proof of 
Proposition^ The cone S is circular with axis gi (dashed line). The cone 7^i 
is bounded by hypeiplanes, and its edges are along the vectors qi,... ,q6. 
The edge qa is the farthest from the axis gi and lies on the surface of S. 


mi 7 ^ and m 2 7 ^ (- 2 , there exists a permutation on 
{ 1 , 2 ,3,4} that maps mi to m 2 and £1 to £ 2 . we see that 
there exists an orthogonal transformation M G "H such that 


and . 

This can be extended to a group action on Q as 




= M 


fmi - 




Ifmi -rrj 

~ V 8 ~ ~ Sm2,^2- 

Thus we see that the group Ti. acts transitively on Q, i.e., 


Q = (Mgj I M G "H} for every i = 1,..., 12. 


and are related to Vi as 

Vi = n 7 ^ 1 , where 7^i = 


The group code Q and the 11 half-spaces Vi can be explicitly 
calculated starting from the geometry of the tetrahedron, and it 
can be verified that 7?,i, and hence Vi, is bounded by exactly 6 
planes arising from 6 of the eleven half-spaces Vi- The region 
72.1 is a convex cone i23i generated from the 6 edges running 
along the vectors qi,..., qe that are the intersections between 
the 6 hyperplanes, i.e., 72i is the infinite cone generated 
from the convex hull of the set {qi,..., qg}. Fig. | 6 ] shows 
an illustration of the geometry considered in this proof (the 
depiction of qi,..., qg is not exact). Since 


T ■ 

mm Ki V = mm -r—r-, 

vGX)i ^ xGT^i X 


(27) 


and since g{x/||x|| is the cosine of the angle between x and 
gi, our problem is to find a vector in 72 i which makes the 
largest angle with gi. The set of points that make a constant 
angle with gi form the surface of an infinite circular cone 
with gi as its axis. Thus dZTl l is equivalent to finding the 
smallest circular cone S, with gi as the axis, that contains the 
conical region 72i. Since 72i is generated by qi,..., qg, S is 
the smallest circular cone that contains the vectors qi,..., qg, 
and has gi as the axis. It follows that S contains on its surface 
the vector q^, from among qi,..., qg, that makes the largest 
angle with gi. Thus, 


j . g{x . g{x 

min gl V = mm 77-^ = min ---^ 
vex>i ^ xG-Ri X xGS X 


This makes Q a group code, and consequently, the regions 
Vi,... , 77 i 2 are congruent to each other ll 22 ll . i.e., for every 
1 < * < j < 12 , there exists an orthogonal transformation 
M G 7£ such that 


The numerical value minjg^i gj giQi / ||qi|| = l/v^ is ob¬ 
tained by a direct computation of the half-spaces T^i,..., T^n, 
and the resulting vectors qi,..., qg arising from the tetrahe¬ 
dral geometry. ■ 


Vj = MVi = {Mv I V G Vi} . 


Since orthogonal transformations conserve inner products and 
since gi G Vi for all i, we have 


min gjV = min g J v for any i 7 ^ 

vGUt vGI>3 ■' 


Thus, to complete the proof it is enough to show that 


min g| V = 
vGKi 



We now restrict ourselves to one particular region Vi and 
find the smallest value of g{v. Note that when v G §^, the 
inner product of v with gi decreases with increasing distance 
IIV — gill . Thus, from (l26T l. Vi is the intersection of with 
the set of all points in that are closer to gi than any other 
gi G G- The region Vi is called a fundamental region of 
the group code G and is bounded by two-dimensional planes 
passing through the origin ll22l . The half-spaces Vi that define 
this fundamental region are 


T’i = {x G I ||x- gill < ||x- g,||} 

= (x G I (gi - g*)'^x > 0} , 


Proposition 2 . If a tetrahedral array is used at the receiver 
and p > 1, then for every receive orientation U, there exists 
fl 2 X 2 submatrix Hgub of the channel matrix H such that 


0 ^ i^suh cos 



where ^sub is the correlation between the two columns of 
Hgub- 


Proof: From Proposition [T] there exist to 7 ^ £ such that 
gJn — •\/l/2. Let Hsub be the submatrix of H formed 


by the to*^ and rows. From 
T] > 1 , we have figuh = 


and the hypothesis that 


<cos(^^). 


The following upper bound on /r* ( 77 ) follows immediately 
from Proposition |2l 


Theorem 2. For a tetrahedral receive array and rj > 1, 

Proof: Let H = [hm,r^ be the 4x2 channel matrix. 
From Proposition m assume without loss of generality that the 
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2x2 submatrix formed from the first two rows has correlation 
Msub < cos (7r/2-\/2?7). Then, 


M = 


< 


2 “b ^2 1^2.2 “b ^3 i^3,2 “b 3^4,2 

h\ 3 / 11,2 + 3 / 1 : ' 

1 


h \ ihz ,2 + /14,1/14,2 


— 2 /^sub + ^ 


/l3,l/l3.2 + /l4,1^4,; 


< - cos 
- 2 






\ 2 V 2 r]) 


+ 


4’ 


where the last inequality follows from Proposition |2] and the 
fact that all hm,n have unit magnitude. ■ 

The upper bound (l + cos (7r/2-\/2)) /2 on /i*(ry) is less 
than 1 for ry > 1. Since d(/i, AX) is a decreasing function of 
fi, we have d(/4*(?7), AX) > d(l, AX). Hence, the geometry 
of the tetrahedral arrangement allows the error probability to 
decay faster than that of rank 1 LoS MIMO channels, and 
provides performance improvement over any planar arrange¬ 
ment Ur = i of antennas, irrespective of the code used at the 
transmitter. Note that this gain of the tetrahedral arrangement 
over planar arrays is not due to larger inter-antenna distances 
dt and dr- 

From (l2Tli . we have d(l, AX) > (||Axi|| — ||Ax 2 ||)^. Us¬ 
ing yi* < 1, we obtain 


d(yt*, AX) > d(l, AX) > (IIAxill - ||Ax2||)' > 0. 


Hence, unlike the planar case, the error probability of a 
tetrahedral receiver is exponential in SNR for any code 

Example 5. Performance of Spatial Multiplexing with Tetra¬ 
hedral Receive Array. Consider the SM scheme signalled 
over rit = 2 antennas using 4-QAM symbols. Let the trans¬ 
mit orientation /3 = 0 be fixed, the inter-terminal distance 
i? = 10 m, A = 4.2 mm, and dt = dr = 0.145 m. 
Then, rj = RX/(2dtdr cos/3) = 1, and from Theorem |2] 
lJ*{v) ^ 0.722. An exhaustive numerical computation over 
all pairs of codewords yields minAx d(0.722, AX) = 0.556. 
Using (l24ll . the pairwise error probability of SM for fixed 
transmit orientation and random receive orientation can be 
upper bounded as 

E(PEP) < PEP* < iexp(-SNRyi*(l)) 

< i exp (—SNR X 0.556). 

On the other hand, as shown in Example |4] for any planar 
receiver array, the error rate is not better than SNR~^. ■ 


B. System design for arbitrary array orientations 

In Section lTV-AI we assumed that rj was fixed, i.e., the trans¬ 
mit orientation /3 and inter-terminal distance R were fixed, and 
we studied the effect of an arbitrary rotation U of the receive 
array on p and error probability. We now design a system that 
allows arbitrary transmit and receive array orientations and a 
range of values Rmin < R < i?max- It is desirable that the 


txi 



Fig. 7. Triangulai* aiTangement of transmit antennas. 


LoS MIMO system guarantees a minimum channel quality i.e., 
M < ftmax, for some /imax < 1- Using (l24l) . for such a system, 

E(PEP) < PEP* < i exp d(yt„,ax, AX)^ . 

Using union bound, the average codeword error rate and bit 
error rate of the system can be upper bounded by 

/ n^SNR , 

— exp I - - — mmd(ytinax, AX) 1 . 

Hence, the coding gain of an arbitrary coding scheme over 
this LoS MIMO system is minAx d(Aimax, AX). 

When the number of transmit antennas n* = 2, by choosing 
/3 = 7r/2, we observe from ( fTOl i that the worst case correlation 
Pmax = 1 irrespective of the array geometry used at the 
receiver. Hence, in order to have ytmax < 1, we need more 
than 2 antennas at the transmitter. 

Suppose the transmitter uses an array of nt > 3 antennas. 
Based on the transmit array orientation, one can choose 2 of 
the Tit antennas for signal transmission so that the angle /3 
corresponding to the chosen pair of antennas is minimum. For 
example, let n* = 3 antennas be placed at the vertices of an 
equilateral triangle with inter-antenna distance dt, as shown in 
Fig. |2] Let tm,n be the unit vector in in the direction of 
the position of transmit antenna m with respect to the position 
of transmit antenna n. Note that the vectors tm,n vary with 
changes in the transmit array orientation. If antennas m and 
n are used for transmission and if u S is the unit vector 
along the direction OO' of transmission, then sin /3 = „ 

(cf. Fig. m where txi and tx 2 correspond to tx^ and tx„, 
respectively). The six vectors in the set 

'T' — {tm,n I tH., Tl — f, 2, 3, m f 

are aiTanged symmetrically in a two-dimensional plane at 
regular angular intervals of 7r/3. Let uy and ua be the 
components of u parallel and perpendicular to the plane of 
T, respectively. Since the vectors in T divide the plane into 
regular conical regions of angular width 7r/3, there exists at 
least one vector tm.n S T such that the angle between tm,n 
and U|| lies in the interval [—7r/6,-|-7r/6], i.e.. 



We can thus upper bound as follows 

|uTt™,„P = + |U|jtm,„P < 0 -f ||U|||pi < i. 
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Fig. 8. The functions fi*, upper bound on fi* and the line /imax = 2/3. 


Thus there exists a tm,n such that |sin(/3)| = < 

1/2, i.e., /3 G [—7r/6,7r/6]. Hence, if the transmit array is an 
equilateral triangle, by appropriately choosing 2 out of the 3 
available antennas for signalling, one can ensure |/3| < tt/6 . 

The upper bound on /r*(ry) of Theorem |2] is not tight and is 
available only for rj > 1. Since this bound can not be used as a 
good estimate of p* (ry) and the analytical computation of the 
exact expression (|2^ of fj,* (p) appears to be difficult, we use 
numerically computed values of p* (p) for system design. The 
function p* (p) and the upper bound of Theorem |2] are shown 
in Fig. m Using the exact function yr*(? 7 ), the requirement on 
channel quality p < ytmax can be translated into a criterion 
r] & [?7min, ??max]- From (l22ll . for fixed dt, dr. A, and \P\ < 
/3max, we have 


ddmin A 

^min — nj 7 ^max — 

Zdfdr 


R„ 


2dtdr C0S/3n 


(28) 


The range ^max] can thus be obtained from ( l28l l. 

Example 6. Suppose we require /Tmax = 2/3 with A = 
4.2 mm. Using a triangular transmit array we have /3niax = 
7r/6. From Fig. [S] the criterion p*{ri) < 2/3 is equivalent to 
??min = Vi = 0-62 and pmax = ?72 = 1-22. If each side of 
the triangular transmit array has length dt = 6 cm, and the 
tetrahedral receive array has dr = 25 cm, then from (l28l) we 
have i?niin = 4.43 m and i?max = 7.75 m. ■ 

The narrow range of [iimin, 7?max] in Example |6] can be 
attributed to the small value of 772 — pi in Fig. |8] This can be 
improved by using a pentagonal transmit array as follows. As 
shown in Fig. |9] with a regular pentagon, the choice of the 
transmit antenna pair can be divided into the following two 
cases: (i) the two antennas are the neighbouring vertices of 
the pentagon with inter-antenna distance equal to the length 
dt of the edge of the regular pentagon, or (ii) the antennas are 
non-neighbouring with inter-antenna distance (l -f dt/2. 
Irrespective of the class from which the antenna pair is 
chosen, it is straightforward to show that |/3| < tt/IO can 
be always guaranteed. While the value of rj for the first case 
is given by (l22l l. in the second case it reduces by a factor of 


(1 -f v^) /2 because of the larger inter-antenna distance. Thus 
the maximum correlation with pentagonal transmit array is 

Fpent(^) = min|p*(p),yr* 

where /r*(p) is given in ( |2^ . From Fig. [8] the value of pmax 
improves from p 2 to pa, thereby widening 7?i„ax]- 

Example 1. As in Example |6] let /tmax = ‘2,t, A = 4.2 mm, 
dt = 6 cm and = 25 cm. With a pentagonal transmit 
array, /Smax = tt/IO, and using the function we 

have pmin = Pi = 0.62 and p^ax = pa = 2. Using (|28l), 
7?min = 4.43 m and i?niax = 12.7 m. ■ 

V. Simulation Results 

We use the system parameters A, dt, dr, i?max and i?min 
from Example |7] We assume that the transmit and receive 
arrays undergo independent uniformly random 3-dimensional 
rotations about their centroids, and the distance R between the 
terminals is uniformly distributed in 7?max]- In all the 

simulations the channel matrix H was synthesized using ([T]l 
and the exact distances {rm,n} between the transmit and 
the receive antennas. We consider the following three coding 
schemes with the transmission rate of 4 bits per channel use: 
(0 the Golden code im using 4-QAM alphabet, (ii) spatial 
multiplexing (SM) ifTSl - lfTTl with 4-QAM, and (Hi) uncoded 
16-QAM transmitted using only one transmit antenna (single¬ 
input multiple-output SIMO). Gray mapping is used at the 
transmitter to map information bits to constellation points, and 
unless otherwise stated, maximum-likelihood (ML) decoding 
is performed at the receiver. While we used pairwise error 
probability for performance analysis in Sections HU HIH and IIVI 
we simulate the bit error rate to compare the average error 
performance. 

A. Error performance with nr = 4 

Fig. [Tol shows the performance of the three schemes with 
two different antenna geometries: (/) uniform linear array 
(ULA) at the transmitter with rtt = 2, and uniform rectan¬ 
gular array (URA) at receive]^ with Ur = 4, (ii) selecting 

^The performance of uniform linear array at receiver is worse than that of 
URA, and hence has been omitted. 


2p 


1 + ^5 
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Fig. 10. Comparison of PentxTetr with ULAxURA. 



Fig. 12. Performance of different tx an'ays with tetrahedral rx an'ay. 



Fig. 11. Coding gain for bit rate of 4 bits per channel use. 


2 antennas from a pentagonal array at the transmitter, and 
using a tetrahedral array at the receiver. The values of dt, 
dr are ideal for the ULAxURA configuration Q at the 
distance R = 2dtdr/X = 7.14 m, which is near the mid-point 
of the interval 7?i„ax]- The performance of the single¬ 

antenna transmission scheme is independent of the receive 
antenna geometry since, from O, all the channel gains 
of the SIMO channel have unit magnitude. Also, Fig. m 
shows the performance of the ideal channel with p, = 0, i.e., 
R = yTv I 2 , which is a pair of parallel AWGN channels each 
carrying a 4-QAM symbol. From Fig. [Tol we see that, with 
ULAxURA, the performance of both SM and the Golden 
code are worse than SIMO at high SNR. Further, since 
minAx d(l 5 AX) = 0 for SM, the error probability decays 
slowly with SNR, confirming our theoretical results. With the 
proposed pentagon x tetrahedron geometry both codes show 
improved performance, close to that of the ideal channel. 

The above error performance is succinctly captured by the 
coding gain minAx d(p, AX) shown in Fig. [TT]as a function 
of p. From Example |7] p < 2/3 for the new antenna geometry. 
From Fig. [TT] we see that the coding gains of SM and the 
Golden code are both equal to 1 for all p < 1/2 and are 
larger than the SIMO coding gain for p < 2/3, which explains 



Fig. 13. En'or probability of spatial multiplexing with triangular transmit 
array when the receive array is (i) three-dimensional, and (ii) rectangular. 
Results are shown for rir = 16 and rir = 64 antennas. 


their superiority to SIMO. On the other hand, the coding 
gain for linear and rectangular arrays is minAx d(l, AX). For 
p = 1, from Fig. [TT] we observe that SIMO has the largest 
coding gain followed by the Golden code and then SM. The 
error performances in Fig. [T0| show this same trend for the 
rectangular array at high SNR. 

Fig.dJjcompares the performance of different transmit array 
geometries when a tetrahedral array is used at the receiver. 
The rit = 2 case (ULA) performs poorly since pmax = 1- 
While the triangular array with the Golden code achieves 
most of the available gain, the pentagonal array has near ideal 
performance. 


B. Error Performance with large number of receive antennas 

The LoS MIMO system analysed in Section|IV|employs the 
tetrahedral receive array - a three-dimensional antenna array 
for Ur = d antennas - to enable smaller error rates than planar 
arrays. The geometry of the receive array is relevant even 
if the number of receiving antennas is large. Theorem [1] 
and Example @1 show that the probability of error of the SM 
scheme is lower bounded up to a constant factor by SNR~^ 
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for any value of n^, if a planar receive array is used. On 
the contrary, from Example |5] the SM scheme can achieve 
exponential rate of decay of error probability if n,. = 4 
antennas are placed at the vertices of a regular tetrahedron. 
It follows that for any > 4, a careful three-dimensional 
arrangement of antennas can ensure that the error rate is 
exponential in SNR. For instance, if the three-dimensional 
arrangement includes a subset of 4 antennas that form a 
tetrahedron, it immediately follows from Example |5] that a 
sub-optimal decoder that bases its decision only on the signals 
received by these 4 antennas achieves exponential error rate. 
Hence, the optimal ML decoder that utilizes all the Ur receive 
antennas achieves an exponential error probability as well. 

Fig. fTSl compares the error performance of SM scheme 
under planar and three-dimensional receive antenna arrays 
when Ur = 16, 64. A triangular array is used at the transmitter, 
4-QAM is chosen as the modulation scheme and ML decoding 
is performed at the receiver. For both values of rir, we consider 
a URA (rectangular arrangement of receive antennas) for the 
planar arrangement of antennas. The three-dimensional array 
is chosen as a set of Ur points on the surface of a sphere so 
that the minimum distance between the points is large. A table 
of such arrangements of points, which are known as spherical 
codes, is available online ll24ll . For fairness, the diameter of 
the sphere is set equal to the width of the rectangular array. 
The coordinates of the points on the sphere were obtained 
from ll24ll . As with previous simulations, we set the values 
of dt, X, Rmax and i?niin as in Example [T] The inter-antenna 
distance dr of the URA is chosen to be 12.5 cm when rir = 16 
and to be 6.25 cm when rir = 64. This is the optimal inter¬ 
antenna distance for the URA when the transmit and receive 
arrays are oriented broadside to each other and the inter¬ 
terminal distance R — l.lAm 0. 

It is evident from Fig.[T3]that array geometry is an important 
design parameter even when rir is large. The error rates of 
rectangular arrays shown in Fig. [13] decay as SNR ^ at high 
SNR. The gain due to the three-dimensional array is about 
7 dB at an error rate of 10“® for both Ur = 16 and 64. 

VI. Conclusion 

We studied the error performance of arbitrary coding 
schemes in 2 x LoS MIMO channels where the commu¬ 
nicating terminals have random orientations. We analyzed the 
effects of some receive array geometries on error probability, 
and showed that, unlike linear, circular and rectangular arrays, 
the error rate with a tetrahedral array decays faster than 
that of a rank 1 channel. Using tetrahedral and polygonal 
arrays, we designed a LoS MIMO system that provides a good 
error performance for all transmit and receive orientations. By 
modelling the R matrix, we derived error probability bounds 
for the case when the number of transmit antennas used for 
signalling is 2. Analysis of the performance when more than 
2 transmit antennas are used is yet to be addressed. 
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